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ABSTRACT 


As  mesoscale  models  increase  in  resolution  there  is  a 
greater  need  to  understand  predictability  on  smaller 
scales.  The  predictability  of  a  model  is  related  to 
forecast  skill.  It  is  possible  that  the  uncertainty  of  one 
scale  of  motion  can  affect  the  other  scales  due  to  the 
nonlinearity  of  the  atmosphere.  Some  suggest  that 
topography  is  one  factor  that  can  lead  to  an  increase  of 
forecast  skill  and  therefore  predictability. 

This  study  examines  the  uncertainty  of  a  mesoscale 
model  and  attempts  to  characterize  the  predictability  of 
the  wind  field.  The  data  collected  is  from  the  summer, 
when  the  synoptic  forcing  is  relatively  benign.  Mesoscale 
Model  5  (MM5)  lagged  forecasts  are  used  to  create  a  three- 
member  ensemble  over  a  12-hour  forecast  cycle.  The 
differences  in  these  forecasts  are  used  to  determine  the 
spread  of  the  wind  field.  Results  show  that  some  mesoscale 
features  have  high  uncertainty  and  others  have  low 
uncertainty,  shedding  light  on  the  potential  predictability 
of  these  features  with  a  mesoscale  model. 

Results  indicate  that  topography  is  a  large  source  of 
uncertainty.  This  is  seen  in  all  data  sets,  contrary  to 
other  studies.  The  ability  of  the  model  to  properly 
forecast  the  diurnal  cycle  also  impacted  substantially  on 
the  character  and  evolution  of  forecast  spread.  The 
persistent  mesoscale  features  were  represented  reasonably 
well,  however  the  detailed  structure  of  these  features  had 
a  fair  amount  of  uncertainty. 
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INTRODUCTION 


I . 


A.  MOTIVATION 

Mesoscale  models  provide  the  user  with  a  more 
realistic  view  of  the  atmosphere,  capturing  details  of 
mesoscale  phenomena  with  higher  resolution  (Doyle  1997  and 
others)  .  However,  this  does  not  necessarily  lead  to  an 
improvement  in  forecast  skill  because  they  are  attempting 
to  model  things  whose  behaviors  and  time  scales  are  not 
fully  understood  or  observed.  In  order  to  have  an  overall 
increase  in  forecast  skill,  the  model  error  for  all 
variables  must  be  reduced  everywhere.  A  study  by  Weygandt 
and  Seaman  (1994)  shows  an  increase  in  model  error  when 
just  horizontal  resolution  is  increased.  The  skill  of  many 
parameters  decreased  suggesting  that  other  factors,  such  as 
model  physics,  vertical  resolution,  and  initial  conditions 
become  more  important  as  grid  resolution  is  enhanced. 

The  impact  of  higher  model  resolution  on  other 
parameters  is  important  because  it  impacts  the  skill  of  the 
model  forecast.  Consequently,  there  becomes  a  greater  need 
to  know  about  predictability  on  small  scales. 

B.  PREDICTABILITY 

Anthes  and  Baumhefner  (1984)  define  predictability  as 
the  upper  limit  to  forecast  skill.  He  suggests  that  an 
inherent  limit  to  the  predictability  of  atmospheric  motions 
exists  because  one  cannot  completely  and  accurately  observe 
the  atmosphere  at  all  times  and  on  all  scales  of  motion. 
As  a  result,  all  forecast  skill  and  predictability  would 
eventually  be  lost  given  enough  time,  due  to  these  inherent 
uncertainties  in  the  initial  conditions. 
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Lorenz  (1982),  another  leader  in  predictability 
studies,  describes  the  bounds  of  predictability.  He 

considered  the  lower  bound  to  be  the  skill  of  current 
forecast  procedures,  while  the  upper  bound  is  based  on  the 
predictability  of  instantaneous  weather  patterns,  or  as 
Lorenz  states,  the  predictability  of  the  instability  of  the 
atmosphere  with  small  amplitude  perturbations.  Lorenz 
(1982)  suggests  that  the  lack  of  predictability  in  the 
behavior  of  the  atmosphere  is  evidence  for  this 

instability.  He  goes  on  to  define  a  predictability  time 
limit.  This  is  the  amount  of  time  between  the  best 
estimate  of  the  atmosphere  based  on  observations  and  an 
estimate  of  its  state  at  a  future  time,  to  the  point  at 
which  the  forecast  looses  all  skill.  After  a  forecast 
reaches  this  limit  it  is  unusable,  and  basically  is  no 
better  than  guessing  (Lorenz  1982).  This  predictability 
limit  is  strongly  dependant  upon  the  accuracy  of  the 
measure  of  the  initial  conditions. 

Figure  1  is  a  graphical  depiction  of  theoretical  error 
growth.  It  compares  the  error  associated  with  climatology, 
persistence  and  the  numerical  model  output.  The  error  in 
the  numerical  model  starts  small  and  then  grows  rapidly 
from  the  initial  time  step  and  is  associated  with  model 
spin  up.  As  the  dynamic  imbalances  are  resolved  the  error 
begins  to  decrease.  There  is  a  period  of  low  error  from 
which  error  begins  to  grow  again  due  to  the  incompleteness 
of  model  physics.  Finally  reaching  a  point  where  the 
forecast  no  longer  has  skill. 

Anthes  (1986)  also  demonstrates  the  growth  of  error  by 
using  the  variance  of  500mb  heights  (Figure  2) .  He 
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compares  the  error  associated  with  a  forecast  of 
climatology  (Ec)  ,  a  persistence  forecast  (Eobs)  ,  a  model 
forecast  (E)  as  seen  in  the  previous  forecast.  He  also 
shows  the  predictability  error  growth  (Ep)  ,  which  represents 
an  estimate  of  the  maximum  obtainable  forecast  skill. 

Lorenz  first  studied  synoptic  scale  predictability  in 
1965.  The  result  of  his  study  demonstrated  that  the  error 
growth  rate  was  dependent  upon  the  synoptic  situation  and 
that  a  7-day  forecast  could  be  possible,  but  a  month  long 
forecast  was  out  of  the  question  (Lorenz  1965)  .  The  result 
of  a  later  study,  by  Lorenz  in  1969,  included  some  small 
scale  effects.  This  study  indicated  that  if  the  large 
scale  was  observed  perfectly,  the  small  scale  uncertainties 
would  induce  error  on  the  large  scale  and  grow  as  if  errors 
existed  in  the  large  scale  initially  (Lorenz  1969)  . 

The  results  of  other  early  predictability  studies 
suggest  that  predictability  would  decrease  as  the 
horizontal  scales  became  smaller.  The  effect  of  synoptic 
pattern  and  weak  versus  strong  instability  was  also  found 
to  alter  the  predictability.  In  summary,  predictability 
had  been  found  to  vary  with  horizontal  scale,  season, 
geographic  location,  and  synoptic  pattern  (Anthes  1986). 
It  was  also  shown  that  when  the  same  set  of  observations 
were  added  to  different  models  they  each  had  different 
error  growths  . 

Anthes  (1986)  confirmed  the  results  of  Lorenz's  1969 
study  also  indicating  that  the  mesoscale  was  less 
predictable  than  the  synoptic  scale.  Anthes  also  suggests 
that  the  uncertainty  of  one  scale  of  motion  would 
contaminate  all  scales  due  to  the  nonlinearity  of  the 
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atmosphere.  Due  to  the  energy  exchanged  among  all  scales 
of  motion  and  the  nonlinear  nature  of  the  atmosphere,  the 
mesoscale  reaches  its  predictability  limit  before  the 
larger  scales.  Tennekes  (1978)  suggested  that  this 
nonlinear  energy  transfer  was  a  result  of  2  and  3- 
dimensional  turbulence,  implying  that  error  on  one  scale 
would  contaminate  the  other  scales  (Anthes  1984) .  Kuypers 
(2000)  found  while  there  was  no  correlation  between 
synoptic  and  mesoscale  error,  the  majority  of  the  error  at 
the  smaller  scales  was  dominated  by  the  error  from  the 
larger  scales. 

Anthes  (1986)  describes  the  predictability  error 
growth  or  representing  the  growth  of  initial  small 
differences  in  a  perfect  model.  He  proposes  that  the 
reason  different  scales  of  motion  were  forecast  with 
different  skill  levels  is  due  to  an  inherent  difference  in 
the  predictability  error  growth  for  each  scale  of  motion. 
This  also  suggests  that  model  error  could  grow  faster  than 
predictability  error,  resulting  in  forecast  skill  less  than 
what  predictability  theory  suggests.  The  estimate  of  this 
model  error  would  be  indicative  of  the  possible 
improvements,  which  could  still  be  realized. 
Contributions  to  model  error  would  be  from 
parameterizat ions ,  numerics,  boundary  conditions  and 
initial  conditions  (Anthes  1986). 

There  is  a  general  impression  that  mesoscale  forecasts 
near  topography  are  more  predictable  (Mass  et  all  2002  and 
others)  .  The  idea  is  that  with  higher  resolution  the 
presence  of  terrain  is  better  represented  helping  to  reduce 
the  error  and  therefore  increase  predictability.  Higher 
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resolution  models  better  simulate  the  topography  and 
therefore  capture  more  of  the  topographically  forced 
features.  A  study  done  by  Weygandt  and  Seaman  (1994) 
showed  a  decrease  of  mean  error  for  topographically  forced 
features,  confirming  that  terrain  increases  predictability. 
The  benefits  to  increasing  model  resolution  for  capturing 
topographically  induced  circulations  were  also  seen  in  a 
study  by  Mass  et  al  (2002) . 

The  other  half  of  this  argument  suggests  that  the 
presence  of  topography  makes  the  mesoscale  less  predictable 
through  an  increase  in  error.  It  was  found  that  the 
largest  wind  speed  errors  in  the  mesoscale  were  near  the 
topography  during  the  case  of  a  landfalling  front  (Nuss  and 
Miller  2001)  .  It  will  also  be  shown  in  this  thesis  that 
the  largest  model  errors  are  associated  with  regions  of 
topography.  The  differences  between  model  topography  and 
reality  can  at  times  be  significant,  contributing  to  error. 

C .  OBJECTIVES 

A  suggested  approach  for  examining  the  predictability 
of  a  mesoscale  model  is  to  determine  the  spread  between  the 
forecasts.  The  idea  being  that  low  spread  events  will  be 
more  predictable  than  those  with  high  spread.  It  has  been 
observed  that  if  the  spread  is  small,  the  skill  of  the 
forecast  tends  to  be  higher  (Steenburgh  2002) . 

Since  the  spread  is  indicative  of  uncertainty  it  is 
used  in  this  study  to  examine  error  growth  and 
distribution.  The  spread  was  determined  using  the 

difference  between  the  various  model  forecasts  and  model 
analyses.  The  amount  of  spread,  or  uncertainty,  and  the 
period  and  location  of  which  it  occurred  are  used  to 
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characterize  the  predictability  of  the  mesoscale  wind 
field . 

The  Mesoscale  Model  version  5  (MM5) ,  run  in  real-time 
by  the  Naval  Postgraduate  School  (NPS)  ,  is  used  to 
construct  lagged  forecast  wind  fields  in  order  to 
investigate  mesoscale  predictability  based  on  a  mesoscale 
model.  The  error  growth  for  each  forecast  cycle  was 
examined  for  the  California  summer,  which  extended  from  01 
May  2002  until  04  October  2002  and  related  to 
predictability  theory. 
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II.  MESOSCALE  MODEL  AND  METHODS 


A.  MODEL  DESCRIPTION 


The  Mesoscale  Model 

version  5 

(MM5 ) 

is  used  for 

all 

model 

runs  and  analyses 

for  this 

thesis  . 

Researchers 

at 

Penn 

State  University 

and  the 

National  Center 

for 

Atmospheric  Research  developed  and  continue  to  support  the 
MM5  (PSU/NCAR  2003)  .  It  is  a  limited-area,  nonhydrostatic, 
terrain-following,  sigma-coordinate  model  designed  to 
simulate  or  predict  mesoscale  and  regional-scale 
atmospheric  circulations.  MM5  can  be  run  on  a  Unix  machine 
or  a  Linux  based  PC.  The  computer  power  needed  to  run  this 
model  increases  as  mesh  size  and  grid  resolution  increase 
(PSU/NCAR  2003)  .  Dr.  Doug  Miller  at  NPS  ran  the  MM5  used 
for  this  thesis. 

The  version  used  in  this  study  is  MM5V2.12,  the  second 
version  of  the  MM5 .  This  version  of  the  model  can  be  run 
hydrostatically  or  non-hydrostat ically  depending  on  what 
the  user  chooses.  The  basis  of  a  hydrostatic  model  is  the 
assumption  of  hydrostatic  equilibrium.  These  models  are 
better  for  global  and  synoptic  features,  since  they  do  not 
account  for  much  of  the  necessary  physics  necessary  to 
model  smaller  features.  Non-hydrostatic  models  basically 
define  a  reference  state  and  the  perturbations  from  that 
state.  These  models  are  preferred  for  mesoscale  modeling 
because  they  account  for  vertical  motions  and  accelerations 
rather  than  inferring  them  from  the  horizontal  convergence 
and  divergence,  i.e.  continuity.  It  becomes  important  to 
use  a  non-hydrostatic  model  when  you  begin  to  talk  about 
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features  that  are  similar  in  height  and  length  scales  and 
on  the  order  of  about  10  km. 

As  with  all  mesoscale  models,  MM5  requires  an  initial 
condition  and  lateral  boundary  condition  for  a  model  run, 
it  is  therefore  coupled  with  a  global  model  or  another 
regional  model.  It  uses  the  other  model's  output  as  a 
first  guess  for  objective  analysis  or  as  the  lateral 

boundary  conditions  (PSU/NCAR  2003) .  In  this  case,  the 
National  Centers  for  Environmental  Prediction  (NCEP) 
Aviation  model  (AVN)  is  the  parent  model  providing  the 
boundary  conditions  in  the  MM5  36-hour  (h)  forecast  as  well 
as  the  analyses  and  forecasts  for  the  basis  of  the  MM5  12-h 
forecast  in  the  case  of  a  cold  start  (Miller  2002) . 

A  model  cold  start  creates  the  analysis  using 
climatology,  or  a  prior  analysis  if  available,  and 
observations  that  have  been  assimilated  into  the  model. 

The  difference  for  a  mesoscale  model  is  that  the  synoptic 
forecast  is  interpolated  down  to  the  mesoscale  and  used 
with  observations.  This  MM5  cold  start  is  generated  from 
the  AVN  analyses  and  forecasts,  as  AVN  is  the  parent 
synoptic  model.  There  were  a  very  small  number  of  cold 
starts  over  the  data  collection  period;  the  majority  were 
warm  starts.  A  warm  start  combines  observations  with  the 
model's  most  recent  forecast,  or  first  guess,  in  order  to 

generate  the  next  forecast.  The  observations  are  used  to 

nudge  the  model  towards  reality.  The  MM5  in  this  study 
uses  the  12-h  forecast  as  the  first  guess  for  the  36-h 
forecast.  Both  types  of  starts  can  induce  errors;  it 
depends  on  the  quality  of  the  parent  field  and  data 
assimilation . 


MM5's  vertical  and  horizontal  resolution  and  domain 
size  are  variable,  allowing  the  user  to  define  them  to 
their  specific  needs  (PSU/NCAR  2003) .  MM5  defines  its 
vertical,  or  sigma,  levels  in  terms  of  pressure.  The  near 
ground  sigma  surfaces  closely  simulate  the  terrain  in  the 
MM5,  while  at  higher-levels  they  tend  to  approximate 
isobaric  surfaces  (PSU/NCAR  2003)  .  The  advantage  to  using 
sigma  over  constant  pressure  or  height  surfaces  is  the  fact 
that  the  sigma  surfaces  do  not  intersect  the  topography  as 
the  other  surfaces  do.  This  allows  the  user  to  easily 
increase  the  vertical  resolution  near  the  surface, 
enhancing  the  representation  in  the  planetary  boundary 
layer  if  desired.  The  number  of  sigma  levels  in  the  MM5 
used  here  is  30  for  all  grids.  Thirteen  of  these  vertical 
levels  are  used  to  describe  the  portion  of  the  atmosphere 
between  the  surface  and  700mb. 

The  MM5  used  in  this  study  is  a  triple-nested  model 
and  therefore  has  three  horizontal  grid  resolutions.  The 
largest  or  the  coarse  resolution  is  108  km  (59  x  59  grid 
points);  next  is  the  fine  resolution  grid  at  36  km  (49  x  61 
grid  points);  and  the  superfine  grid  at  12  km  (91  x  127 
grid  points)  is  the  smallest.  These  grids  are  nested  over 
the  west  coast  of  the  United  States  and  California  as  seen 
in  Figure  3 . 

B.  VISUAL  AND  GEMPAK 

VISUAL  is  a  FORTRAN-coded  diagnostic  and  display 
program  that  uses  NCAR  graphics  utility  routines  to  look  at 
meteorological  information  (Nuss  and  Drake  1995).  This 
program  was  used  to  examine  the  statistical  data  set  at 
various  pressure  levels  as  well  as  vertical  cross  sections. 
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VISUAL  also  allows  for  generated  plots  to  be  printed  for 
publication . 

The  General  Meteorological  Package,  or  GEMPAK,  is  a 
suite  of  application  programs  for  the  analysis,  display, 
and  diagnosis  of  meteorological  data  (UCAR  1998) . 

C .  PROCEDURE 

MM5  raw  forecast  data,  in  the  form  of  GEMPAK  files, 
were  collected  beginning  01  May  2002  through  04  October 
2002,  basically  encompassing  the  California  dry  season. 
Each  model  run  went  out  to  a  36-h  forecast,  in  3-h 
increments.  In  order  to  ensure  a  three  member  ensemble, 
lagged  forecasts  were  used.  Lagged  forecasting  used  the 
24-h,  12-h  and  initial  forecasts,  which  were  all  valid  for 

the  same  time  period,  based  on  the  assumption  that  there 
would  not  be  a  large  difference  between  them.  In  a 
relatively  unchanging  synoptic  regime,  such  as  the 
California  summer,  this  was  a  valid  assumption.  Due  to  the 
use  of  lagged  forecasts,  only  the  periods  up  to  the  12-h 
forecast  could  be  used.  Only  one  or  two  members  would  be 
present  in  the  ensemble  if  we  had  tried  to  use  forecasts 
beyond  the  12-h  forecast.  Statistics  were  calculated 
during  the  data  processing.  The  mean  wind  speed  and  spread 
(difference  in  wind  speed  from  the  forecasts  and  associated 
analysis)  were  calculated.  The  spread  becomes  important 
because  we  believe  this  is  indicative  of  model  error. 

The  GEMPAK  files  were  first  interpolated  from  the 
108km  and  36km  domains  onto  the  12km  domain.  No  artificial 
effects  resulted  from  the  interpolation.  This  step  of 
processing  the  data  also  allowed  for  it  to  be  further 
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processed  for  statistics  and  plotted  using  the  VISUAL 
program . 

Next  was  an  attempt  to  get  a  basic  plan  of  how  to 
categorize  the  synoptic  regimes  to  create  smaller  groups. 
The  data  was  plotted  in  VISUAL  using  500mb  heights  and 
850mb  winds  for  each  day  at  0000UTC  and  1200UTC  (hereafter 
OOz  and  12z  model  analysis)  .  It  was  determined  to  divide 
the  season  into  four  categories  based  on  the  850mb  wind 
pattern  to  investigate  the  impact  of  different  synoptic 
flow  patterns.  The  four  categories  were:  along  coast  flow 
(coast  parallel),  onshore  flow,  offshore  flow,  and  weak 
flow.  The  flow  over  the  region  from  Cape  Mendocino  to  Big 
Sur  was  used  to  determine  the  analysis  category.  The  data 
set  was  also  examined  as  one  large  group. 

Each  day  was  then  processed  to  determine  the  mean  and 
spread  of  the  wind  field  over  the  12  km  domain.  The  mean 

was  the  mean  wind  speed  and  the  spread  was  the  maximum 

difference  in  the  wind  speed  between  the  three  lagged 
forecasts.  The  program  took  three  12-h  forecast  periods 
with  five  forecasts  for  each  period  (3  hour  forecast 
intervals)  .  Table  1  is  an  example  of  how  the  lagged 
forecasts  were  analyzed.  In  this  example  the  program  used 
the  24-h  forecast  from  the  OOz  analysis,  the  12-h  forecast 
from  the  12z  analysis  on  the  30th  and  the  OOz  analysis  on 

the  1st.  The  program  used  those  times  and  then  went  out  to 

the  next  forecast,  in  three-hour  intervals. 
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Step 

1 

30/00z  f 2 4 

30/12z  f 12 

Ol/OOz  fOO 

Step 

2 

f  2  7 

f  15 

f  03 

Step 

3 

f  30 

f  18 

f  0  6 

Step 

4 

f  33 

f  2 1 

f  0  9 

Step 

5 

f  36 

f  2  4 

f  12 

Table  1:  Example  of  how  the  statistical  analysis  program, 
calc-stats,  processed  the  forecast  data.  In  the  form  of: 
30/00z  f24  =  30 (day) /  OOz  (model  analysis)  f24  (24  hour 

forecast) . 

The  next  step  was  to  average  all  of  the  mean  data  and 
all  of  the  spread  data  for  each  time  period  and  forecast. 
This  was  done  for  each  flow  pattern  (along,  onshore, 
offshore  and  weak)  as  well  as  for  the  entire  data  set.  An 
input  file,  which  consisted  of  each  date  and  analysis  time, 
was  first  created  for  the  flow  patterns  and  then  processed 
for  the  average.  These  data  sets  were  not  broken  down  by 
model  start  time,  so  OOz  and  12z  were  averaged  together. 
In  order  to  create  a  manageable  input  file  for  the  whole 
season  each  month  was  averaged  by  analysis  time  (i.e.  each 
month  had  two  files  a  OOz  and  a  12z) .  The  months  were  then 
averaged  together  into  larger  input  files  representing  the 
OOz  and  12z  for  the  summer  season.  The  processing  was  done 
in  this  manner  due  to  the  amount  of  data  files  and  does  not 
affect  the  results  since  it  is  a  simple  averaging  of  data. 

This  data  was  further  analyzed  with  plots  created  in 
VISUAL.  The  levels  850mb  and  500mb  were  plotted  to 

describe  the  synoptic  patterns  over  the  entire  domain  at 
those  levels.  All  of  the  plots  created  were  from  the 
statistical  data,  the  mean  wind  speed  and  its  spread.  The 
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plots  created  were  for  the  seasonal  data  as  well  as  the 
four  flow  categories.  Four  West-East  cross  sections  were 
chosen  as  illustrated  in  Figure  4.  The  cross  sections  were 
chosen  for  their  differences  in  terrain.  The  vertical 
cross  sections  were  plotted  for  each  model  start  time,  OOz 
and  12z,  as  well  as  for  the  00  to  12  hour  forecast  of  each 
run  (forecasts  are  every  3  hours)  .  The  vertical  cross 
sections  for  the  four  flow  patterns  were  also  plotted  by 
forecast  time.  Since  these  forecast  times  were  averaged 
together  in  the  previous  step  they  contain  both  OOz  and  12z 
model  start  times  and  therefore  have  only  one  set  of 
forecasts  (not  two  like  the  summer  case) .  Horizontal  plots 
were  also  created  in  a  similar  fashion  for  each  of  the  data 
sets . 

The  final  step  used  to  analyze  the  data  sets  was  to 
construct  an  average  of  the  spread  over  the  entire  domain 
for  each  sigma  level.  The  output  of  this  program  was  by 
level,  from  the  surface  to  500mb  (21  levels)  ,  and  the  mean 
value  of  the  spread  field  for  each  forecast  (to  12  hours) 
of  the  OOz  and  12z  model  start  times.  This  program  was 
initially  run  using  the  entire  12  km  domain.  However,  we 
found  this  to  be  strongly  influenced  by  the  large  area  over 
the  ocean,  where  the  spread  decreased  through  the  12-h 
period  independent  of  the  synoptic  situation.  These  level 
averages  were  not  representative  of  the  regions  near 
topography  where  the  spread  changed  significantly  with  each 
forecast  as  seen  in  the  VISUAL  plots.  The  increased  spread 
near  the  topography  was  essentially  masked  by  the  low 
spread  values  over  the  water  in  the  average.  In  an  attempt 
to  better  illustrate  the  over  land  performance,  the 
averaging  domain  was  reduced  to  limit  the  amount  of  data 
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over  the  water  and  focus  on  the  coastal  and  inland 
topographical  effects.  This  new  domain  considered  the 
eastern  two-thirds  of  the  12  km  domain  and  was  therefore  no 
longer  dominated  by  the  trends  over  the  ocean. 
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III.  GENERAL  DISCUSSION 


A.  SUMMER  CLIMATOLOGY 

The  climate  of  California  is  quite  varied;  along  the 
coast  it  is  mild  and  tends  to  be  cooler  in  the  northern  and 
central  parts  of  the  state,  while  the  southeastern  region 
is  hot  and  dry.  California's  summer  season,  also  the  dry 
season,  begins  in  the  late  spring  and  continues  until 
October  when  the  wet  season  begins.  The  primary  synoptic 
feature  during  this  time  period  is  the  East  Pacific 
(EASTPAC)  High.  This  dominates  the  region  with  large-scale 
subsidence  and  produces  northwesterly  surface  winds  over 
the  eastern  Pacific  Ocean.  These  northwesterlies  drive  the 
Ekman  transport  along  the  west  coast  of  North  America 
forcing  the  surface  waters  offshore.  This  allows  the  cold 
deep  water  to  rise  to  the  surface  resulting  in  very  cool 
ocean  temperatures  near  the  coast.  The  strong  subsidence 
and  cool  ocean  temperatures  result  in  a  strong  low-level 
inversion  and  a  well  developed  boundary  layer  through  out 
the  region. 

A  thermally  induced  surface  low  over  the  desert 
Southwest  and  the  interior  region  of  California  is  also 
present  during  the  summer  months  as  a  result  of  significant 
daytime  heating.  This  thermal  trough  helps  establish  a 
strong  cross-coast  pressure  gradient  and  strong 
northwesterly  flow. 

B.  OVERVIEW  OF  DYNAMICS 

There  are  several  mesoscale  events  that  occur  during 
the  summer  season.  The  prominent  features  were  seen  in  the 
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analysis  of  this  thesis.  This  section  describes  the  basic 
dynamics  of  these  features. 

1 .  Coastal  Jet 

The  California  coastal  jet  is  present  a  majority  of 
days  during  the  summer  months.  The  jet  flows  along  the 
coast  towards  the  south  with  the  strongest  winds  near  the 
coast  and  weaker  winds  offshore.  Its  maximum  is  neither  at 
the  surface  nor  up  against  the  coastal  mountains,  but 
rather  slightly  elevated  and  away  from  the  coastal  terrain 
(where  the  maximum  pressure  and  temperature  gradients  are 
located) .  This  is  due  to  surface  friction,  which  retards 
the  wind  closer  to  the  surface.  The  jet  is  initiated  by 
the  cross-coast  pressure  gradient  established  by  the 
synoptic  regime  and  is  forced  by  the  coastal  topography. 
Wind  speeds  in  the  jet  core  can  be  greater  than  40  knots 
with  little  diurnal  variation. 

There  is  a  strong  baroclinic  structure  at  the  coast 
and  a  well-mixed  marine  boundary  layer  offshore.  The  depth 
of  the  inversion  increases  offshore  due  to  the  weakening  of 
the  synoptic  subsidence  and  the  warmer  sea  surface 
temperatures.  The  isentropes  slope  with  the  inversion 
leading  to  the  coastal  baroclinic  structure  and,  through 
thermal  wind  relation,  leads  to  a  low-level  wind  maximum. 
It  has  been  observed  that  the  stronger  the  isobaric  slope 
the  stronger  the  jet  can  become  (Nuss  2002) . 

Coastal  topography  also  plays  a  role  in  jet  dynamics. 
The  northwesterly  flow  is  typically  blocked  by  the  terrain 
due  to  the  presence  of  a  strong  low-level  inversion  and  is 
turned  down  gradient  parallel  to  the  coast.  The  flow  is 
strongest  near  the  coast  and  weakens  offshore  as  it 
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approaches  the  Rossby  radius  of  the  mountains.  This  will 
occur  even  if  there  is  no  slope  to  the  inversion.  This  is 
also  consistent  with  the  structure  implied  by  the  thermal 
gradient  (Nuss  2002)  . 

2 .  Mountain  and  Valley  Breezes 

The  mountain  and  valley  winds  are  diurnally  forced. 
They  are  generally  favored  when  there  is  weak  synoptic  flow 
and  a  weak  pressure  gradient.  Formation  of  these  breezes 
depends  upon  the  contrast  in  surface  temperatures,  the 
difference  between  daytime  heating  and  nighttime  cooling, 
the  orientation  of  the  mountain  slope  (heating  is  strongest 
on  the  southern  and  eastern  facing  slopes) ,  and  the 
direction  of  the  synoptic  flow.  It  is  also  important  to 
have  mostly  clear  skies  in  order  for  the  strongest  heating 
and  cooling  to  occur.  (COMET  2003) 

a .  Mountain  Breeze 

The  mountain  breeze  consists  of  an  up-slope  and 
down-slope  wind.  When  the  sun  rises  it  begins  to  heat  the 
slopes  of  the  mountain  and  the  air  above  causing  the  warmer 
air  to  rise  up  the  slope.  The  result  of  this  up-slope  wind 
is  subsidence  in  the  valley.  By  the  afternoon  the  up-slope 
winds  are  at  their  peak.  As  the  solar  heating  decreases 
and  radiative  cooling  begins  the  winds  near  the  surface 
begin  to  reverse  while  those  higher  in  the  boundary  layer 
still  remain  up-slope  winds.  As  the  mountains  cool  the 
whole  system  reverses  due  to  the  development  of  cooler, 
denser  air.  This  action  results  in  the  down-slope  winds. 
The  air  in  the  valley  rises  in  response  to  the  down-slope 
winds  (Ahrens  1994)  . 
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b.  Valley  Breeze 

The  valley  breeze  occurs  along  with  the  mountain 
breeze.  The  up-slope  winds  begin  shortly  after  sunrise  and 
after  sufficient  heating  has  occurred  in  the  valley.  Later 
in  the  morning  the  winds  become  up-valley  winds.  The  wider 
and  deeper  the  valley,  the  longer  it  takes  for  the  winds  to 
shift  direction.  The  up-valley  breeze  lasts  until  after 
sunset.  Once  adequate  cooling  has  taken  place  the  down- 
valley  breeze  takes  over.  This  lasts  until  after  sunrise 
when  the  cycle  begins  again.  The  valley  breezes  occupy  the 
lower  10  to  30%  of  the  total  valley  depth.  The  average 
up/down-valley  breeze  speed  is  10  m/s,  but  can  be  stronger 
depending  on  the  strength  and  depth  of  the  inversion. 
These  winds  tend  to  accelerate  as  they  travel  through  the 
valley.  As  they  exit  the  valley  they  gradually  slow  as  the 
flow  spreads  out  (COMET  2003) . 

3.  Thermal  Effects 

Thermal  circulations  are  very  similar  to  sea  breeze 
circulations,  but  involve  the  difference  of  heating  over 
land  instead  of  the  land  and  water  temperature  differences. 
If  there  is  no  variation  in  temperature  (or  pressure)  in 
the  horizontal  over  the  land  there  will  not  be  a 
circulation.  In  this  situation  isobaric  surfaces  are  flat. 
When  the  atmosphere  is  heated  in  one  area  more  than 
another,  or  cooled  more  in  one  area  versus  another,  the 
isobaric  surfaces  become  sloped.  They  are  close  together 
over  the  cooler  region  and  spread  apart  over  the  warmer 
region.  This  slope  leads  to  a  pressure  gradient  aloft, 
with  the  air  moving  from  high  to  low  pressure.  There  is  an 
associated  vertical  circulation  with  the  horizontal  one. 
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caused  by  warm  air  rising  and  cold  air  sinking.  The 
surface  lows  and  highs  created  through  this  process  are 
referred  to  as  thermal  highs  and  lows. 

4 .  Mountain  Waves 

Mountain  waves  are  created  when  the  wind  blows  over  a 
ridge;  the  air  parcels  are  displaced  vertically  and,  if  the 
air  is  stably  stratified,  will  descend  and  may  oscillate 
about  their  equilibrium  level.  The  result  is  a  gravity 
wave,  referred  to  as  a  mountain  or  lee  wave.  For  ranges 
that  are  less  than  100  km  in  width  the  perturbations  seen 
in  the  wind  fields  are  primarily  in  the  vertical,  where  as 
ranges  with  widths  greater  than  100  km  the  perturbations 
are  predominately  in  the  horizontal  (Durran  1986) . 

The  initial  perturbation  of  the  air  parcels  triggered 
by  the  mountains  was  seen  most  predominately  in  the  Central 
California  cross-section  due  to  interaction  with  the 
Sierra-Nevada  Mountains.  This  range  in  eastern  California 
is  about  690  km  (430  miles)  long  and  60-110  km  (40-70 
miles)  wide,  with  many  of  its  peaks  reaching  4270  m  (14,000 
ft)  or  higher.  The  perturbation  of  the  air  parcels  in  this 
case  is  in  the  vertical  due  to  the  width  of  the  range. 

The  stability  of  the  atmosphere  plays  a  role  in  the 
formation  and  appearance  of  the  mountain  waves.  These 
waves  are  buoyancy  driven  oscillations  and  therefore  need  a 
stably  stratified  environment  in  order  to  support 
propagation.  If  a  weaker,  stable  layer  is  above  a 
stronger,  stable  layer  the  wave  amplitude  will  decay  with 
height.  The  waves  could  also  be  trapped  in  this  case.  The 
period  of  these  oscillations,  or  the  Brunt-Vaisala 
frequency  (N2)  ,  and  the  frequency  imposed  by  the  terrain  as 
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the  air  flows  over  are  also  important  in  determining  the 
vertical  propagation,  the  amplitude,  and  trapping  of  these 
waves . 
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IV.  RESULTS 


A.  OVERVIEW  OF  RESULTS 

The  structure  and  characteristics  of  mesoscale 
uncertainty  are  illustrated  in  four  cross  sections  taken 
through  the  12  km  domain  (Figure  4)  .  The  northern  most 
cross  section  proceeds  through  Cape  Mendocino  and  the 
northern  parts  of  the  Sierra-Nevada  mountain  range;  the 
next  cross  section  goes  through  San  Francisco  and  the 
central  region  of  the  Sierra-Nevada  range;  the  third  cross 
section  runs  near  Point  Conception  toward  Lake  Havasu;  and 
the  final,  most  southern  cross  section  extends  through  San 
Clemente  Island,  just  north  of  San  Diego  into  Arizona. 
Each  cross  section  was  chosen  due  to  the  differences  in 
topography.  All  cross  sections  start  offshore,  in  order  to 
capture  the  majority  of  the  coastal  jet  feature,  and  end 
just  east  of  the  California  border.  The  persistent 
mesoscale  features  seen  in  this  study  will  be  examined 
further.  These  features  include  the  coastal  jet,  mountain 
and  valley  circulations,  topographic  effects,  and  the 
structure  of  the  model  atmosphere  with  height. 

The  expected  trend  of  error  in  wind  speed  to  grow  over 
time  is  not  observed  in  this  study  (Figure  1)  .  Although 
the  spread  did  increase  between  the  analysis  and  the  3-h 
forecast  of  each  forecast  cycle,  the  spread  did  not 
significantly  decrease  and  then  slowly  increase  as  might  be 
expected  as  the  model  adjusts  to  the  initial  conditions. 
Instead,  the  spread  of  the  wind  speed  tended  to  remain 
rather  large  through  the  12-h  forecast  cycle  as  observed  in 
the  level-average  tables,  (Tables  2-7)  .  This  is 
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particularly  true  at  lower  levels,  while  in  all  cases  the 
spread  increased  in  the  upper  levels  throughout  the  12-h 

forecast  cycle.  There  is  also  a  difference  in  the  OOz  and 
12z  forecast  cycles.  The  OOz  cycle  has  significantly  more 
spread  through  the  12-h  forecast  than  the  12z  cycle  does, 

at  all  levels. 

A  diurnal  cycle  is  also  observed  in  the  model.  The 

OOz  forecast  series  have  more  uncertainty  associated  with 
all  of  the  features  than  the  12z  cycle  does.  This  is  seen 
in  all  of  the  forecast  hours.  There  is  an  increase  in 

uncertainty  in  the  3-h  forecast  for  both,  however  the  OOz 
forecast  cycle  continues  to  have  large  spread  values 
through  the  remaining  part  of  the  forecast.  The  12z 
forecast  cycle  does  not  continue  to  have  large  spread 
values  after  the  3-h  forecast. 

The  primary  relationship  between  spread  and  mean  wind 
speed  found  in  this  study  is  that  as  mean  wind  speed 
increases,  the  spread  also  increases.  This  is  seen  in  all 
cross  sections,  as  an  increase  in  spread  with  height,  and 
is  associated  with  the  larger  mean  wind  speeds  due  to  the 
stronger  dynamics  aloft.  In  addition,  the  mesoscale 

features  with  higher  wind  speeds,  such  as  the  coastal  jet, 
also  have  higher  spread  as  observed  in  the  cross  sections. 

B .  COASTAL  JET 

The  mean  wind  fields  in  MM5  clearly  depict  the  coastal 
jet  throughout  the  summer  season  (Figures  5-27)  .  The  jet 
extends  from  Northern  California  where  it  is  most 
prevalent,  to  Southern  California.  Using  spread  as  an 
indicator  of  uncertainty,  the  intensity  and  position  of  the 
jet  core  exhibited  considerable  variation.  The  mean  jet 
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tended  to  be  lower  in  the  atmosphere  and  further  to  the 
west  than  the  position  of  the  spread  maximum.  In  addition, 
the  spread  was  larger  at  higher  mean  wind  speeds. 

The  cross  sections  demonstrate  the  vertical  structure 
of  the  coastal  jet.  The  strongest  winds  are  closest  to  the 
coast  and  weaken  offshore.  Although  the  jet  typically 
slopes  up  with  the  inversion  from  the  coast  as  it  moves 
offshore  this  slope  is  barely  evident  in  the  cross  sections 
represented  in  Figures  5  and  6.  The  solid  lines  in  these 
figures  represent  the  total  mean  wind  speed,  and  the  dashed 
lines  represent  the  spread.  This  convention  is  used 
throughout  the  rest  of  the  vertical  cross  sections  and 
horizontal  plots.  The  coastal  jet  has  little  diurnal 
variation;  a  comparison  of  Figures  5  and  6  show  a  slightly 
stronger  jet  at  OOz  f03  than  at  the  12z  f03  forecast. 
Figures  5  and  6  show  that  the  jet  core  is  strongest  just 
after  the  maximum  heating  and  is  weaker  and  further  from 
the  coast  during  the  cool  part  of  the  diurnal  cycle.  The 
spread  however,  remains  close  to  the  coast  during  these 
periods.  In  fact,  the  position  of  the  spread  maxima  for 
all  of  the  cross  sections  and  times  remained  close  to  the 
coast.  It  did  not  change  position  with  the  slight  movement 
of  the  jet  core. 

The  magnitude  of  the  jet  is  the  strongest  in  the 
northern  most  cross  section,  which  has  a  range  of  mean  wind 
speed  values  of  ll-13m/s  and  spread  values  of  3.5-4.0m/s. 
The  central  cross  section  was  slightly  lower  with  mean  wind 
speed  values  of  10-12m/s  and  3.0-4.0m/s  for  the  spread.  As 
expected,  the  numbers  continued  to  decrease  as  one  moved 
further  south,  since  the  coastal  jet  tends  to  weaken 
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towards  the  south.  The  values  at  the  Point  Conception 
cross  section  were  9-10m/s  for  the  mean  wind  speed  and  2.5- 
3 . 5m/ s  for  the  spread.  The  weakest  of  all  is  the  Southern 
California  cross  section.  It  has  a  range  of  mean  wind 
speed  values  of  8-9m/s,  with  the  spread  having  a  range  of 
2.5-3. 5m/ s . 

Examining  the  time  variation  of  the  jet  structure  and 
uncertainty  reveals  that  the  coastal  jet  is  stronger  in  the 
OOz  forecast  cycle  than  in  the  12z  cycle.  The  OOz  forecast 
cycle  in  the  Northern  cross  section  is  strongest  at  the  3-h 
forecast  with  a  mean  wind  speed  of  13  m/s  and  spread  of 
4 . 5m/ s  (Figure  5),  while  the  12z  is  at  its  weakest  at  the 
6-hour  forecast  with  a  mean  wind  speed  of  11  m/s  and  spread 
of  3 . 5m/s .  A  reason  for  this  is  the  OOz  run  begins  at  the 
warmest  part  of  the  day,  4pm  local  time.  The  jet  is  well 
represented  in  the  observations  that  are  assimilated  into 
the  model  for  that  run.  The  jet  in  the  OOz  forecast  cycle 
increases  at  the  3-h  forecast  from  12m/ s  to  13m/ s  with  a 
spread  increase  of  0.5m/s,  from  4.0  to  4.5  m/s .  The 
coastal  jet  weakens  after  the  OOz  f03  forecast  for  the 
remainder  of  the  forecast  period.  The  only  exception  is 
the  12-h  forecast  in  the  Point  Conception  cross-section 
where  there  is  an  increase  in  jet  speed  and  spread. 

This  pattern  of  jet  intensity  and  spread  for  the  OOz 
and  12z  forecast  cycles  varied  to  some  extent  with 
location.  The  magnitude  of  jet  speed  and  spread  for  the  3- 
h  forecast  of  the  OOz  forecast  cycle,  is  not  the  strongest 
in  all  of  the  cross  sections.  The  cross  section  near  Point 
Conception  has  its  strongest  OOz  jet  at  the  3-h  forecast. 
The  other  cross  sections  have  their  maxima  at  the  6-h 


forecast  for  the  OOz  series.  The  weakest  jet  speed  in  the 
12z  run  is  consistently  the  6-h  forecast,  with  the 
exception  of  the  Southern  California  cross  section  where 
the  weakest  is  at  the  initial  time. 

The  850mb  horizontal  plots  show  the  varying  horizontal 
structure  of  the  coastal  jet  in  the  different  synoptic  flow 
regimes  as  well  as  the  small  changes  in  the  summer  case  as 
a  result  of  diurnal  effects.  These  plots  are  not  ideal,  in 
that  they  do  not  go  through  the  maximum  of  the  coastal  jet, 
but  instead  pass  right  above  it.  However,  the  850mb  winds 
can  affect  the  intensity  of  the  coastal  jet  and  are 
generally  indicative  of  its  horizontal  location  (Nuss 
2002)  . 

The  OOz  and  12z  forecasts  show  that  the  coastal  jet  is 
strongest  in  Northern  California,  weakens  around  the 
Central  Coast  and  continues  to  weaken  as  one  looks  south 
(Figure  7)  .  The  along  coast  flow  shows  the  coastal  jet 
extending  strongly  into  the  Southern  California  region,  as 
noted  by  the  increased  wind  speeds  along  the  coast  (Figure 
8)  .  In  all  of  the  other  flow  cases  the  coastal  jet  begins 
to  relax  further  to  the  north.  The  offshore  flow.  Figure 
9,  is  very  similar  to  the  OOz  and  12z  jet.  The  jet  is 
slightly  stronger  offshore  (8  vs.  lOm/s)  as  it  extends  down 
the  coast  along  115W,  but  is  still  weaker  in  Southern 
California.  The  onshore  flow  does  not  have  a  strong  jet 
with  mean  winds  less  than  5m/ s  even  in  the  northern  part  of 
the  state  (Figure  10)  .  A  weak  increase  in  mean  wind  speed 
as  one  moves  offshore  is  apparent  and  the  only  indication 
of  a  coastal  jet.  The  maximum  mean  wind  speed  of  the  jet 
at  this  level  is  7m/s,  where  the  previous  cases  were  9m/ s 
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or  greater.  The  weak  flow  has  a  very  prominent  jet  in 
Northern  California,  with  mean  wind  speeds  over  8m/s,  but 
the  flow  rapidly  weakens  to  less  than  3m/ s  just  north  of 

San  Francisco  Bay  (Figure  11)  .  The  weak  flow  case  has  this 
weakening  of  the  flow  along  the  coast  further  north  than 

any  other  case. 

Variations  in  the  coastal  jet  and  its  associated 
spread  under  the  four  different  synoptic  flows  are  also 
evident  in  the  cross  sections.  For  example,  in  the  along 
coast  case,  the  coastal  jet  is  stronger  in  Southern 
California  than  in  any  other  flow  regime  (Figure  12)  .  The 
spread  in  this  region  is  no  larger  than  in  the  summer 
average,  despite  the  increase  in  mean  wind  speeds.  The 

location  of  the  spread  also  remains  near  the  coast.  The 
constant  spread  values  for  the  along  coast  flow  regime 
compared  to  the  summer  average,  suggest  that  MM5  is  less 
sensitive  to  synoptic  variability  in  this  flow  regime. 

The  weak  flow  coastal  jet  had  very  large  spread 

values.  It  demonstrates  large  uncertainty,  which  is  seen 
in  the  spread  as  mean  wind  speed  increases.  It  also  shows 
the  inability  of  the  model  to  forecast  the  position  and 
intensity  in  a  weak  flow  situation.  In  this  flow  regime, 
the  Southern  California  region  had  the  weakest  jet.  The 
uncertainty  here  is  more  from  variation  in  intensity  than 
in  position,  since  the  spread  contours  match  up  relatively 
well  with  those  of  the  mean  wind  speed  (Figure  13) . 

The  along  coast  flow  cases  do  not  have  the  strongest 
coastal  jet  in  all  cross  sections  as  might  be  expected. 
The  Cape  Mendocino  section  has  a  mean  wind  speed  value  of 
13m/s  for  all  forecast  times  with  spread  values  ranging 
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from  2 . 5m/ s  to  3.25m/s.  Nevertheless,  these  values  are  not 
the  largest.  The  weak  flow  Cape  Mendocino  case  is  the 
strongest  and  also  has  the  largest  uncertainty  associated 
with  it.  It  has  a  constant  mean  wind  speed  of  14m/s  over 
all  forecast  periods  and  a  spread  ranging  from  4 . 5m/s  to 
5m/ s .  The  offshore  flow  Cape  Mendocino  case  is  also  larger 
than  the  along  coast,  with  a  mean  wind  speed  of  14m/s  for 
all  forecast  times  and  a  spread  of  3.0-3.25m/s. 

C.  MOUNTAIN  AND  VALLEY  CIRCULATIONS 

Mountain  and  valley  breezes  are  seen  in  all  of  the 
cross  sections  to  varying  degrees.  The  northern  cross 
sections  have  both,  as  does  the  Point  Conception  profile, 
while  the  Southern  California  cross  section  has  a  mountain 
breeze . 

In  the  northern  most  cross  section,  a  valley  breeze  is 
noticed  between  the  coastal  mountains  and  the  Sierra- 
Nevadas .  This  feature  exhibits  a  lot  of  variability 
(Figure  5)  .  During  the  OOz  forecast  cycle  the  mean  wind 
speed  of  this  breeze  is  between  5  and  8m/s,  with  the 
strongest  times  at  the  OOz  f03  and  OOz  f06  forecasts.  The 
spread  throughout  the  OOz  forecasts  varies  as  well.  The 
range  was  3.0-4.5m/s,  yet  the  largest  spread  does  not 
always  correspond  to  the  larger  mean  wind  speeds  here  as 
observed  in  other  features.  The  large  spread  associated 
with  this  feature  most  likely  occurs  during  the  transition 
from  an  up-valley  to  a  down-valley  breeze.  This  switch 
transpires  after  diurnal  heating  and  cooling  begins.  The 
12z  cycle  starts  off  with  a  mean  wind  speed  of  6m/ s  and  a 
spread  of  3.5m/s,  as  the  forecast  time  increases  the  mean 
and  spread  of  the  wind  field  decreased  to  4m/s  and  2.5m/s, 
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respectively.  MM5  has  difficulties  with  the  cooling  part 
of  the  diurnal  cycle,  which  adds  to  the  uncertainty  in  the 
forecast  skill  of  this  feature.  This  breeze  is  one  of  the 
weaker  features,  with  low  mean  wind  speeds,  but  the 
associated  spread  is  large.  The  greatest  spread  also 
occurs  at  the  time  when  the  model  is  staying  warm  instead 
of  cooling  in  the  OOz  forecast  cycle.  These  same 
characterist ics  are  observed  for  all  of  the  different 
synoptic  flow  cases  as  well. 

The  12z  f03  forecast  of  the  San  Francisco  Bay  cross 
section  shows  increased  uncertainty  in  the  region  near  the 
steep  slope  of  the  topography,  which  is  associated  with  a 
mountain  breeze  (Figure  14)  .  This  feature  is  not  very 
strong,  but  is  highly  variable.  Through  the  rest  of  the 
12z  forecast  the  feature  is  weak  but  predictable  since  the 
uncertainty  is  very  low.  Unlike  the  coastal  range,  which 
has  a  high  amount  of  uncertainty  throughout  both  forecast 
periods,  the  increased  spread  is  only  on  the  windward  side 
at  the  one  forecast  time.  MM5  may  not  be  fully  resolving 
the  steepness  of  the  windward  slope  in  this  region.  This 
uncertainty  is  present  during  a  forecast  from  the  cool  part 
of  the  cycle  (12z  f03)  ,  which  MM5  handles  well.  MM5  may 
over  heat  this  western  facing  slope  too  early  causing  a 
difference  in  the  transition  from  down-slope  to  up-slope 
winds  as  well  as  in  mean  wind  speed. 

This  mountain  breeze  is  also  seen  in  all  of  the  cases 
for  the  different  synoptic  flows.  As  with  the  summer  12z 
case  (Figure  14),  all  of  the  flows  similarly  show  a  spread 
maximum  at  the  3-h  forecast  that  decreases  significantly  in 
the  later  forecasts.  All  of  the  cases  have  mean  wind 
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speeds  of  4m/s  in  that  region.  The  weak  and  offshore  flows 
have  a  spread  of  3.5m/s,  which  is  greater  than  the  3 .  Om/s 
spread  of  the  along  coast  and  onshore  flows. 

The  lee  side  of  the  coastal  mountains  in  the  Southern 
California  cross  section  has  a  region  of  high  spread  in  all 
forecasts  of  the  summer  case  and  in  all  of  the  various 
synoptic  flow  cases  (Figure  15)  .  The  largest  spread  is  in 
the  initial  and  3-h  forecasts  for  all  cases.  There  is  a 
trend  for  the  maximum  spread  to  weaken  beginning  with  the 
6-h  forecast  in  all  cases  as  well.  This  area  of  large 
uncertainty  is  probably  due  to  the  model  not  correctly 
portraying  the  mountain  breeze. 

D .  TOPOGRAPHIC  EFFECTS 

The  characteristics  in  the  four  different  synoptic 
flow  cases  are  surprisingly  similar  to  each  other  and  to 
the  summer  average.  This  is  not  expected  since  the  regimes 
are  all  very  different;  along  coast,  offshore,  onshore  and 
weak  flow.  It  was  anticipated  that  the  mesoscale  features 
noted  would  be  different  for  each  flow  pattern  due  to  their 
potential  differences  in  interaction  with  topography. 

The  large  region  of  high  spread  seen  on  all  of  the 
Southern  California  cross  sections  on  the  lee  side  of  the 
coastal  range  is  a  result  of  MM5  incorrectly  forecasting  a 
mountain  breeze  (Figure  15)  .  This  feature  is  present  at 
all  time  steps  and  cases  suggesting  that  MM5  over  forecasts 
the  up  and  down-slope  winds  in  this  region.  The  mean  winds 
are  rather  light  throughout  the  season,  but  the  spread  here 
is  exceptionally  high,  especially  during  the  OOz  forecasts, 
when  the  model  should  have  been  cooling  the  atmosphere. 
The  uncertainty  during  the  warming  cycle  is  large  here  as 
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well,  but  less  than  during  the  cooling  cycle  (>5m/s  vs. 
3.5m/s) .  The  model's  diurnal  problem  most  likely  adds  to 
the  uncertainty  seen  in  the  OOz  run.  Another  possible 
explanation  for  this  large  spread  is  that  MM5  is  not 
showing  the  across  mountain  flow  properly.  This  flow  is 
due  to  the  cool  temperatures  over  the  water  and  the  warm 
temperatures  over  land  to  the  east  of  this  range.  It  is 
also  possible  that  the  model  is  not  correctly  forecasting 
the  low  level  stratification  allowing  a  mountain  breeze  to 
take  place  when  in  fact  it  is  being  prevented  by  the 
stratification . 

The  fluctuations  seen  in  the  wind  fields  over  the 
mountains  are  presumably  due  to  flow  interaction  with  the 
topography  and  could  be  representative  of  mountain  wave 
formation  in  the  model  (Figure  16)  .  Durran  (1986)  showed 
that  for  steeper  slopes  the  wave  amplitude  is  larger,  as  is 
the  tendency  for  vertical  propagation.  This  is  seen  in  the 
different  cross  sections  as  well.  In  general  the  spread  in 
the  upper  levels  is  higher,  but  in  the  regions  of  these 
waves  the  local  spread  is  no  larger  than  the  environment. 
The  spread  fields  are  perturbed  as  well,  so  in  fact,  the 
spread  over  the  mountains  is  actually  lower  than  the  spread 
of  the  surrounding  atmosphere  at  a  given  level. 

In  the  central  cross  section  (Figure  16)  the  terrain 
has  an  effect  on  the  mean  wind  field  that  extends  above 
500mb.  This  cross  section  is  the  only  one  with  significant 
topography;  the  eastern  portion  is  where  it  crosses  the 
Sierra-Nevada  Mountains.  The  mountains  cause  a 

perturbation  in  the  mean  and  the  spread  over  the  peak. 
This  perturbation  is  a  result  of  the  air  flowing  over  the 
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mountains  and  probably  the  mark  of  mountain  wave  formation. 
The  northern  and  southern  cross  sections  also  show 
perturbations  in  the  spread  and  mean  wind  fields,  though 
not  to  the  extent  of  the  central  section.  The  effect  of 
the  coastal  topography  in  those  sections  reaches  only  to 
the  650mb  level. 

All  of  the  four  synoptic  flow  cases  for  the  Central 

California  cross  section  show  perturbations  in  the  mean 
wind  field  and  spread  above  500mb  over  the  Sierra-Nevada 
Mountains  as  in  the  summer  plots  (Figure  16)  .  Similar  to 
the  summer  case,  these  waves  are  the  response  of  the  flow 
as  it  runs  into  the  topography  and  are  presumably  an 
indicator  of  mountain  wave  formation.  The  along  coast  and 
onshore  cases  show  the  largest  amplitude  of  these  waves 
(Figure  17)  .  The  offshore  case  demonstrates  smaller 

amplitude  waves  which  decrease  quickly  with  height  (Figure 
18) .  The  weak  flow  case  has  the  perturbations  as  well,  but 
the  amplitudes  are  not  as  great  (not  shown) . 

The  offshore  flow  cross  sections  are  very  similar  to 
the  summer  average  cross  sections  as  far  as  seeing 

perturbations  associated  with  mountain  wave  formation.  In 
all  of  these  cross  sections  the  perturbations  stay  below 
the  650mb  level.  This  is  true  for  the  weak  flow  cross 
sections  as  well.  The  along  coast  case  does  as  well, 

except  for  the  Point  Conception  region  where  very  small 
amplitude  perturbations  do  occur  near  the  500mb  level.  The 
amplitudes  of  the  waves  decrease  with  height  and  the  lower 
level  disturbances  do  extend  downstream  (Figure  19)  .  The 
weak  flow  Point  Conception  region  is  an  exception  since  the 
waves  in  the  mean  wind  and  spread  propagate  above  500mb. 
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The  mean  flow  in  this  region  was  the  slowest  with  speeds 
between  4  and  6m/ s . 

The  onshore  flow  case  in  the  Southern  California  cross 
section  shows  perturbations  from  the  coastal  topography 
extending  above  500mb  (Figure  20)  .  The  amplitude  of  the 
perturbations  decreases  with  height  and  continues 
downstream.  This  is  different  than  all  of  the  remaining 
Southern  California  cross  sections  where  the  perturbations 
due  to  topography  remain  below  the  600mb  level  and  may  be 
due  to  the  stronger  cross  mountain  flow  in  this  case. 

E .  DIURNAL  VARIATION 

The  portion  of  the  diurnal  cycle  sampled  over  the  12 
hours  impacts  the  evolution  of  uncertainty.  In  the  OOz 
forecast  cycle,  the  spread  grows  significantly  until  the 
12-h  forecast.  Conversely,  the  spread  in  the  12z  cycle 
does  not  increase  as  much  until  the  9  and  12-h  forecasts 
(lpm  and  4pm  local  time)  .  OOz  occurs  at  the  maximum  of 
heating  for  the  day,  4pm  local  time.  It  is  suggested  that 
the  model  tries  to  keep  the  atmosphere  warm  longer,  having 
difficulties  with  the  cooling  portion  of  the  diurnal  cycle. 
12z  is  at  the  cool  part  of  the  diurnal  cycle  and  the 
increase  in  spread  at  the  9  and  12-h  forecasts  demonstrates 
the  ability  of  MM5  to  better  represent  the  warming  portion 
of  the  diurnal  cycle.  This  trend  reveals  the  weakness  of 
MM5  to  cool  the  atmosphere  in  the  OOz  forecast  cycle.  The 
four  synoptic  flow  cases  are  not  separated  by  model  run  and 
therefore  the  diurnal  effects  are  averaged  together.  These 
effects  are  not  observed  in  this  manner  for  those  cases. 

The  Point  Conception  cross  section  illustrates  this 
effect  of  model  start  time  with  respect  to  the  diurnal 
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cycle  (Figures  21  and  22)  .  The  terrain  in  this  cross 
section  is  such  that  some  surfaces  begin  to  heat  before 
others,  which  results  in  a  thermal  circulation.  The  OOz 
run  shows  that  the  model  maintains  a  stronger  circulation 
since  it  does  not  cool  down  very  well.  This  is  manifested 
by  an  increase  in  spread  from  the  fOO  to  f06  forecasts. 
The  12z  run  handles  the  warming  part  of  the  diurnal  cycle 
well.  MM5  picked  up  on  this  circulation.  The  error  is  low 
and  wind  speeds  are  fairly  strong.  Once  the  cooling 
begins,  the  temperature  gradient  relaxes,  the  mean  wind 
speeds  decrease  and  the  spread  increases.  The  higher 
spread  at  this  time  is  due  not  only  to  the  slow  cooling  of 
the  model  but  also  due  to  the  inherent  lack  of 
predictability  of  light  and  variable  winds. 

Diurnal  effects  do  not  considerably  affect  the  coastal 
jet.  However,  MM5  showed  that  the  jet  is  stronger  during 
the  warmer  parts  of  the  day.  This  suggests  that  the  model 
is  able  to  represent  the  cross-coast  thermal  gradient  and 
the  inversion  relatively  well.  As  noted  in  the  earlier 
section  on  the  coastal  jet,  the  jet  movement  in  towards  the 
shore  and  away  from  it,  as  well  as  intensity  differences, 
are  due  to  the  model  representation  of  the  diurnal  cycle. 

F.  STRUCTURE  WITH  HEIGHT 

The  spread  of  the  wind  speed  is  found  to  increase  with 
height,  which  is  a  result  of  the  increase  in  wind  speed 
with  height  due  to  stronger  dynamics  aloft.  The  average 
spread  was  3 . Om/ s  or  greater  at  heights  above  600mb  in  all 
cross  sections.  This  was  also  seen  on  the  level  average 
tables  for  all  (Tables  2-7). 
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In  all  of  the  vertical  cross  sections  it  is  observed 
that  as  one  moves  up  in  the  atmosphere,  the  spread 
increases.  At  the  higher  levels  the  synoptic  scale  should 
have  more  of  an  influence,  implying  that  the  winds  at 
higher  levels  should  be  forecast  with  less  uncertainty. 
Since  the  opposite  is  seen  it  suggests  that  the  mesoscale 
12km  grid  is  responding  to  error  growth  from  the  synoptic 
scale.  This  is  also  supported  by  the  fact  that  the  500mb 

spread  is  greatest  at  all  forecast  times  and  for  all  cases. 

The  spread  at  500mb  also  increases  as  forecast  time 
increases,  which  is  consistent  with  synoptic  scale  error 
growth.  The  synoptic  scale  should  dominate  at  the  longer 
forecast  times,  especially  at  the  higher  levels.  This  is 
also  seen  in  the  level  average  tables.  These  trends  will 
most  likely  increase  if  the  forecasts  are  run  for  longer 
periods  of  time. 

Scale  interaction  is  noted  in  the  spread  values  in  the 
plots  and  in  the  level  averages  (Tables  2-7)  .  There  is 
evidence  of  large  spread  in  the  upper  levels  suggesting 

that  the  synoptic  scale  error  is  beginning  to  affect  the 

model  level  by  level.  The  growth  is  seen  in  level  21 
(500mb)  and  propagates  down  with  increasing  time.  This 
trend  will  probably  continue  out  past  the  12-h  forecast,  if 
this  study  were  able  to  do  so. 

In  the  lower  levels  spread  decreases  after  the  initial 
forecast  period  increase,  yet  had  we  been  able  to  see 
longer  forecasts  we  probably  would  have  observed  an 
increase  at  these  levels,  as  the  synoptic  scale  error  would 
begin  to  dominate  the  forecast.  At  the  same  time  the  error 
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from  the  surface  infiltrates  the  levels  above  it  and  can  be 
seen  slowly  increasing  with  forecast  time. 

The  500mb  plots  show  no  interaction  with  topography, 
though  the  perturbations  due  to  the  vertically  propagating 
waves  over  the  Sierra-Nevada  range  could  be  seen  in  some 
plots  (Figure  23)  .  The  mean  wind  speeds  and  spread  were 
largest  in  these  plots.  The  OOz  forecast  cycle  continues 
to  increase  in  spread  throughout  the  12-h  period.  The  12z 
increases  for  the  first  two  forecasts,  decreases  at  12z  f06 
and  12z  f09  finally  increasing  at  the  last  forecast.  As 

expected,  there  is  no  evidence  of  a  coastal  jet  in  any  of 
the  500mb  plots.  The  mean  winds  in  the  flow  cases  do  show 
evidence  of  the  direction  of  the  synoptic  forcing  in  each 
group  (Figures  24-27) .  The  along  shore  flow  experienced 
the  strongest  forcing  from  the  northwest,  since  this  is 
where  the  strongest  mean  wind  speeds  are  located  (Figure 
24) .  The  offshore  flow  cases  have  the  significantly 
stronger  mean  wind  speeds  to  the  north,  as  did  the  onshore 
flow  case  (Figure  25  and  26)  .  However  the  onshore  flow 
case  has  stronger  mean  wind  speeds.  The  strongest  mean 
wind  speeds  in  the  weak  flow  case  occur  in  the  northeastern 
region  of  the  domain,  suggesting  this  is  where  the  dominant 
synoptic  forcing  is  located  (Figure  27) . 

The  spread  at  the  low  levels,  where  the  flow  interacts 
with  terrain,  is  much  larger  than  it  is  at  the  higher 
levels.  This  is  seen  in  all  of  the  plots  as  definite 
features  near  the  topography.  The  spread  in  those 

mesoscale  regions  frequently  is  larger  than  2.5m/s.  The 
spread  at  the  upper  levels  is  greater  than  2 . 5m/s  above 
600mb  over  a  broad  region  for  all  cases  and  forecast 
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periods.  This  large  low  level  spread  is  not  seen  in  the 
statistics  (Tables  2-7)  .  This  is  due  to  the  averaging  of 
the  entire  level.  These  large  uncertainties  are  in  such  a 
small  part  of  the  domain  that  they  basically  get  masked  in 
the  average  by  the  other  regions  in  the  domain  with  low 
spread  values.  This  indicates  that  low=level  uncertainty 
is  tied  to  mesoscale  features,  while  upper=level 
uncertainty  is  more  tied  to  synoptic  scale  features. 
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VI.  DISCUSSION  AND  CONCLUSIONS 


A.  DISCUSSION 

It  was  noted  that  a  common  characterist ic  of  mesoscale 
models,  the  exaggeration  of  thermal  circulations  with 
topography,  also  exists  in  MM5 .  MM5  seems  to  force  a 
mountain  breeze  even  though  it  may  not  actually  occur. 
This  could  be  due  to  the  fact  that  MM5  is  having  trouble 
accurately  forecasting  the  low  level  stratification  which 
would  inhibit  such  a  strong  mountain  breeze  as  illustrated 
in  the  Southern  California  cross  sections.  It  is  also 
possible  that  MM5  has  difficulty  with  steep  mountain 
slopes.  Recall  that  MM5  is  a  terrain  following  model  so  it 
sees  the  winds  on  the  slope  of  the  mountain  as  horizontal 
winds.  It  then  tries  to  force  the  circulation  from  the 
high-pressure  valley  to  the  low-pressure  mountaintop 
despite  the  environmental  characteristics,  which  may  impede 
the  actual  formation  of  this  breeze. 

1 .  Topography 

The  thought  that  topography  imparts  greater 
predictability  is  not  shown  in  this  study.  Whether  it  is 
the  physical  representation  of  the  terrain  in  the  model  or 
the  circulations  that  result  from  its  presence,  it  is 
illustrated  here  that  the  largest  source  of  uncertainty  is 
the  topography.  Figure  28  depicts  MM5's  version  of 
California's  terrain. 

The  terrain  in  the  model  is  a  lot  smoother  than  the 
real  terrain.  MM5  also  sees  mountains  as  steps  rather  than 
a  steep  slope.  There  are  differences  in  the  model 
topography  as  well.  A  barrier  may  or  may  not  exist  to  the 
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extent  that  it  does  in  reality,  preventing  flows  from  being 
represented  correctly  and  increasing  error.  The  studies 
that  demonstrated  improved  skill  due  to  topography  may  have 
had  a  better  terrain  representation  or  improved  model 
physics . 

The  results  of  this  study  suggest  that  near  and  around 
the  topography  are  where  the  lowest  skill  is  and  the 
greatest  reduction  of  the  predictability  of  this  model. 
This  result  is  also  supported  by  the  low  and  rather 
constant  spread  values  over  the  Pacific  Ocean.  There  is  a 
higher  level  of  predictability  here  than  over  land  and  less 
mesoscale  structure  in  general.  The  valleys  also  had  large 
changes  in  spread.  MM5  may  not  be  forecasting  the  valley 
breezes  properly,  which  is  also  due  to  differences  in  the 
representation  in  terrain. 

2 .  Diurnal  Cycle 

The  results  of  this  study  show  that  the  model  does  not 
accurately  portray  the  diurnal  cycle,  and  that  is  another 
large  source  of  uncertainty.  The  observed  patterns  of 
error  growth  and  decay  support  this.  The  12z  run  follows 
the  trend  that  one  would  expect  and  is  reflected  in  Anthes' 
error  curve  (Figure  2)  .  The  spread  growth  in  the  OOz 
forecast  cycle  does  not  resemble  the  error  curve  in  Figure 
1  or  Figure  2 .  The  uncertainty  seen  here  was  in  part  due 
to  the  model's  inability  to  fully  capture  the  diurnal 
cycle.  MM5  did  not  cool  as  fast  as  the  atmosphere  keeping 
the  model  warmer  longer.  This  lead  to  large  spread  values 
through  the  forecast  cycle. 
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B.  CONCLUSIONS 

Error  growth  is  affected  by  many  different  variables. 
In  this  study  we  only  explored  the  wind  fields  our  attempt 
to  provide  insight  into  the  predictability  of  a  mesoscale 
model.  On  the  whole,  the  basic  characteristics  of  the 
summer  season  and  the  flow  cases  are  very  similar.  This 
therefore  suggests  that  the  overall  uncertainty  of  primary 
mesoscale  features  in  this  model  is  not  that  large.  This 
is  significant  when  trying  to  determine  the  predictability 
of  features  in  a  model.  The  consistency  in  capturing  the 
dominant  mesoscale  features  is  encouraging,  although  the 
predictability  of  the  detailed  structure  can  be  quite  large 

The  standard  error  growth  with  respect  to  time,  as 
seen  in  Figures  1  and  2,  was  not  observed  in  this  study. 
The  initial  uncertainties  were  largest  in  the  first  time 
step.  Upon  further  examination  of  the  level  averages  we 
saw  that  there  is  no  one  point  in  time  that  continually  has 
the  minimum  spread  values.  The  location  of  the  minimum 
spread  changes  with  location  and  case.  This  suggests  that 
either  no  model  adjustment  to  the  initial  conditions 
occurred  or  it  occurred  prior  to  the  3-h  forecast.  More 
likely  is  that  the  large  spread  at  the  analysis  time  masks 
any  mesoscale  error  growth.  This  limits  our  ability  to 
assess  actual  predictability  error  growth,  but  clearly 
highlights  the  uncertainty  in  mesoscale  structure  and  the 
inability  of  the  model  to  correct  for  this  uncertainty. 

C.  SOURCES  OF  ERROR 

There  were  sources  of  error  in  the  way  that  we  chose 
to  approach  this  study.  The  lagged  forecasts  while 
indicative  of  synoptic  scale  variation  were  not  able  to 
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adequately  represent  small  mesoscale  perturbations  with  the 
12-h  time  lag.  A  mesoscale  ensemble  with  small  mesoscale 
perturbations  would  be  a  better  approach.  However,  the 
inability  to  consistently  represent  mesoscale  details  at 
the  initial  times  in  the  model,  found  in  this  study,  is 
problematic  for  any  ensemble  approach. 

The  region  that  was  used  to  determine  the  flow 
categories  left  out  Southern  California.  Hence,  those 
cross  sections  were  not  always  in  a  flow  defined  by  the 
category  they  were  placed  in.  This  also  affected  the 
inland  regions.  The  flow  affecting  the  valley  breeze  could 
have  been  oriented  in  such  a  way  as  to  result  in  position 
or  intensity  changes.  These  flows  were  not  grouped 
together.  The  flows  were  also  not  indicative  of  the  flow 
experienced  by  the  inland  mountains,  such  as  the  Sierra- 
Nevada  range . 

These  errors  could  be  removed  by  averaging  multiple 
summer  seasons  and  would  potentially  provide  better  insight 
to  the  predictability  of  this  model. 

D .  FURTHER  STUDY 

In  order  to  make  this  study  complete,  several  other 
aspects  should  be  considered  which  provide  opportunities 
for  further  research.  Several  are  listed  here. 

1 .  Research  the  Other  Seasons 

Other  seasons  need  to  be  examined  in  order  to  see  what 
happens  during  the  transition  seasons  and  in  the  winter. 
Other  mesoscale  features  will  arise;  there  will  be  more 
synoptic  scale  forcing  events  (summer  was  quite 
uneventful);  and  error  trends  will  vary  due  to  those 
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differences.  Other  weather  patterns  could  provide 

different  results. 

2 .  Study  Numerous  Summer  Seasons 

Collect  data  over  several  summers.  A  more  or  less 
predictable  feature  could  mark  each  summer  skewing  the 
results.  Averaging  these  summers  would  lead  to  a  better 
idea  of  error  and  predictability. 

3.  Research  Different  Model  Parameters 

It  is  also  important  to  study  parameters  other  than 
wind.  This  should  be  done  for  this  data  set  as  well  as  for 
other  seasons.  Model  error  growth  and  predictability 
patterns  found  here  based  on  winds  may  not  hold  for  all 
parameters . 

4 .  Topography 

Terrain  played  a  large  part  in  the  error  of  this 
study,  due  to  the  model  terrain  being  smoother  than 
reality.  Improving  model  topography  would  provide 

different,  possibly  more  accurate,  error  results. 

5 .  Compare  with  Observations 

Observations  should  be  compared  to  the  data  set  to  see 
if  the  features  noted  are  still  valid.  The  perturbations 
over  the  mountains  in  the  wind  field  are  present  through 
out  all  the  cases.  Observations  would  show  if  this 
actually  occurs  as  constantly  as  the  model  suggests. 

6 .  Additional  Statistical  Techniques 

It  would  also  be  useful  to  examine  the  data  set  using 
other  statistics.  This  may  shed  more  light  on  the  subject 
providing  different  results.  The  very  basics  were  used 
here,  mean  and  spread,  and  still  provided  some  useful 
results . 
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APPENDIX  A.  TABLES 


Level 

Pressure 

(mb) 

OOz  fOO 

OOz  f03 

OOz  f06 

OOz  f09 

OOz f 1 2 

1 

22m  AGL 

2.214 

2.002 

1.994 

1.942 

1.908 

2 

1010 

2.389 

2.422 

2.331 

2.218 

2.174 

3 

1000 

2.544 

2.637 

2.522 

2.375 

2.314 

4 

990 

2.623 

2.854 

2.734 

2.558 

2.478 

5 

980 

2.672 

3.000 

2.911 

2.700 

2.599 

6 

970 

2.672 

3.060 

3.001 

2.778 

2.654 

7 

960 

2.639 

3.063 

3.027 

2.806 

2.665 

8 

950 

2.614 

3.061 

3.041 

2.825 

2.671 

9 

925 

2.615 

2.962 

2.962 

2.784 

2.645 

10 

900 

2.496 

2.881 

2.889 

2.724 

2.614 

11 

875 

2.473 

2.839 

2.862 

2.706 

2.611 

12 

850 

2.523 

2.863 

2.913 

2.748 

2.651 

13 

825 

2.480 

2.899 

2.994 

2.856 

2.731 

14 

800 

2.523 

2.916 

3.015 

2.905 

2.773 

15 

775 

2.570 

2.885 

2.978 

2.877 

2.772 

16 

750 

2.609 

2.859 

2.948 

2.863 

2.781 

17 

700 

2.777 

2.916 

2.973 

2.946 

2.905 

18 

650 

2.947 

3.089 

3.128 

3.140 

3.136 

19 

600 

3.242 

3.341 

3.365 

3.371 

3.397 

20 

550 

3.480 

3.570 

3.582 

3.632 

3.697 

21 

500 

3.728 

3.765 

3.837 

3.959 

4.047 

Table  2.  Level  Average:  Summer  OOz 

Average  spread  for  the  12km  domain  at  each  level  and 
forecast  period. 

AGL :  Above  Ground  Level 
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Level 

1 2z  fOO 

1 2z  f03 

1 2z  f06 

1 2z  f09 

12z  f  1 2 

1 

2.406 

2.152 

1.878 

1.773 

1.802 

2 

2.424 

2.292 

2.147 

2.020 

1.989 

3 

2.632 

2.452 

2.230 

2.097 

2.070 

4 

2.729 

2.606 

2.297 

2.169 

2.167 

5 

2.870 

2.715 

2.324 

2.193 

2.225 

6 

2.962 

2.772 

2.325 

2.178 

2.232 

7 

2.988 

2.801 

2.311 

2.134 

2.202 

8 

3.006 

2.813 

2.305 

2.099 

2.175 

9 

3.013 

2.801 

2.327 

2.068 

2.134 

10 

2.895 

2.793 

2.363 

2.087 

2.127 

11 

2.875 

2.799 

2.394 

2.126 

2.144 

12 

2.942 

2.796 

2.393 

2.154 

2.177 

13 

2.943 

2.810 

2.445 

2.223 

2.246 

14 

2.932 

2.810 

2.527 

2.310 

2.320 

15 

2.879 

2.798 

2.604 

2.394 

2.380 

16 

2.828 

2.803 

2.678 

2.477 

2.439 

17 

2.919 

2.885 

2.817 

2.655 

2.605 

18 

3.017 

3.023 

2.983 

2.895 

2.889 

19 

3.201 

3.246 

3.217 

3.169 

3.187 

20 

3.436 

3.484 

3.475 

3.411 

3.443 

21 

3.720 

3.736 

3.713 

3.649 

3.693 

Table  3.  Level  Average:  Summer  12z. 
Defined  as  in  Table  2. 
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Level 

fOO 

f03 

f06 

f09 

f  1 2 

1 

2.130 

1.861 

1.756 

1.668 

1.644 

2 

2.116 

2.068 

1.994 

1.818 

1.721 

3 

2.281 

2.260 

2.131 

1.929 

1.827 

4 

2.378 

2.445 

2.269 

2.057 

1.962 

5 

2.479 

2.590 

2.387 

2.164 

2.076 

6 

2.560 

2.692 

2.460 

2.231 

2.149 

7 

2.581 

2.752 

2.498 

2.254 

2.166 

8 

2.577 

2.775 

2.512 

2.260 

2.168 

9 

2.609 

2.740 

2.488 

2.244 

2.165 

10 

2.506 

2.694 

2.473 

2.237 

2.153 

11 

2.502 

2.665 

2.491 

2.257 

2.172 

12 

2.614 

2.659 

2.512 

2.297 

2.220 

13 

2.554 

2.670 

2.569 

2.391 

2.301 

14 

2.540 

2.681 

2.604 

2.447 

2.369 

15 

2.534 

2.648 

2.602 

2.461 

2.405 

16 

2.518 

2.614 

2.596 

2.472 

2.424 

17 

2.604 

2.612 

2.599 

2.522 

2.510 

18 

2.658 

2.693 

2.669 

2.651 

2.692 

19 

2.808 

2.856 

2.835 

2.820 

2.864 

20 

2.967 

3.007 

3.005 

3.004 

3.059 

21 

3.140 

3.142 

3.209 

3.247 

3.348 

Table  4.  Level  Average:  Along  Coast  Flow. 
Defined  as  in  Table  2. 
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Level 

fOO 

f03 

f06 

f09 

f  1 2 

1 

2.126 

1.926 

1.768 

1.674 

1.675 

2 

2.149 

2.106 

1.942 

1.823 

1.855 

3 

2.316 

2.297 

2.082 

1.949 

1.960 

4 

2.388 

2.481 

2.233 

2.088 

2.086 

5 

2.482 

2.597 

2.344 

2.184 

2.172 

6 

2.518 

2.645 

2.398 

2.229 

2.207 

7 

2.520 

2.657 

2.411 

2.234 

2.210 

8 

2.521 

2.674 

2.424 

2.238 

2.214 

9 

2.551 

2.657 

2.437 

2.243 

2.226 

10 

2.444 

2.628 

2.440 

2.245 

2.233 

11 

2.449 

2.618 

2.441 

2.255 

2.234 

12 

2.529 

2.628 

2.459 

2.272 

2.254 

13 

2.521 

2.652 

2.510 

2.329 

2.306 

14 

2.525 

2.652 

2.551 

2.379 

2.348 

15 

2.510 

2.636 

2.574 

2.412 

2.371 

16 

2.498 

2.639 

2.605 

2.461 

2.409 

17 

2.595 

2.709 

2.706 

2.585 

2.549 

18 

2.699 

2.810 

2.825 

2.747 

2.761 

19 

2.906 

2.995 

2.982 

2.949 

2.994 

20 

3.140 

3.224 

3.178 

3.178 

3.249 

21 

3.407 

3.467 

3.411 

3.452 

3.505 

Table  5.  Level  Average:  Offshore  Flow. 
Defined  as  in  Table  2. 
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Level 

fOO 

f03 

f06 

f09 

f  1 2 

1 

2.061 

1.872 

1.771 

1.723 

1.713 

2 

2.021 

1.932 

1.941 

1.942 

1.834 

3 

2.265 

2.079 

2.061 

2.066 

1.968 

4 

2.310 

2.254 

2.181 

2.209 

2.133 

5 

2.402 

2.375 

2.267 

2.292 

2.263 

6 

2.455 

2.423 

2.297 

2.302 

2.301 

7 

2.450 

2.429 

2.288 

2.259 

2.273 

8 

2.460 

2.432 

2.291 

2.229 

2.246 

9 

2.509 

2.421 

2.257 

2.150 

2.179 

10 

2.420 

2.428 

2.253 

2.102 

2.130 

11 

2.415 

2.473 

2.291 

2.135 

2.148 

12 

2.456 

2.554 

2.380 

2.218 

2.229 

13 

2.470 

2.607 

2.488 

2.335 

2.340 

14 

2.517 

2.635 

2.565 

2.435 

2.417 

15 

2.523 

2.619 

2.599 

2.488 

2.447 

16 

2.504 

2.597 

2.623 

2.529 

2.470 

17 

2.646 

2.644 

2.675 

2.650 

2.582 

18 

2.736 

2.780 

2.812 

2.842 

2.793 

19 

2.881 

2.963 

2.992 

3.032 

3.019 

20 

3.012 

3.087 

3.107 

3.168 

3.217 

21 

3.102 

3.128 

3.193 

3.296 

3.368 

Table  6.  Level  Average:  Onshore  Flow. 
Defined  as  in  Table  2. 
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Level 

fOO 

f03 

f06 

f09 

f  1 2 

1 

2.149 

1.973 

1.759 

1.672 

1.671 

2 

2.091 

2.067 

1.873 

1.800 

1.825 

3 

2.257 

2.270 

2.017 

1.892 

1.916 

4 

2.427 

2.492 

2.173 

2.005 

2.040 

5 

2.546 

2.657 

2.281 

2.070 

2.089 

6 

2.588 

2.717 

2.325 

2.090 

2.069 

7 

2.595 

2.720 

2.331 

2.090 

2.056 

8 

2.638 

2.723 

2.347 

2.102 

2.064 

9 

2.710 

2.675 

2.387 

2.123 

2.054 

10 

2.571 

2.658 

2.394 

2.128 

2.064 

11 

2.501 

2.643 

2.372 

2.120 

2.077 

12 

2.478 

2.614 

2.351 

2.123 

2.088 

13 

2.421 

2.583 

2.399 

2.197 

2.140 

14 

2.423 

2.537 

2.442 

2.259 

2.166 

15 

2.390 

2.496 

2.468 

2.265 

2.162 

16 

2.362 

2.477 

2.508 

2.290 

2.189 

17 

2.413 

2.522 

2.589 

2.477 

2.325 

18 

2.483 

2.612 

2.710 

2.682 

2.571 

19 

2.639 

2.759 

2.883 

2.892 

2.870 

20 

2.749 

2.873 

3.004 

3.036 

3.068 

21 

2.885 

2.964 

3.048 

3.120 

3.187 

Table  7.  Level  Average:  Weak  Flow. 
Defined  as  in  Table  2. 
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APPENDIX  B.  FIGURES 


The  following  pages  of  figures  are  grouped  together  in 
this  appendix  in  order  to  help  in  the  reading  of  this  work. 
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ERROR  VARIANCE  (m2) 


DAYS 


Figure  2.  Growth  of  Error  Variance.  (From:  Anthes 
1986) 


Figure  3.  Model  nested  grid  and  domain  sizes.  (From: 
Miller  2003) 
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Figure  4.  Cross  Section  Plot. 
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Figure  12.  Southern  California  Cross  Section: 
Along  Coast  Flow  fOO. 

Mean  wind  speed  in  m/s  (solid)  and  spread  in  m/s 
(dashed) . 


Figure  13.  Southern  California  Cross  Section: 
Weak  Flow  f06. 

Mean  wind  speed  in  m/s  (solid)  and  spread  in  m/s 
(dashed) . 


Figure  14.  Central  California  Cross  Section:  12z  f03 
Mean  wind  speed  in  m/s  (solid)  and  spread  in  m/s 
(dashed) . 
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Figure  15.  Southern  California  Cross  Section:  OOz  f03 
Mean  wind  speed  in  m/s  (solid)  and  spread  in  m/s 
(dashed) . 


Figure  16.  Central  California  Cross  Section:  12z  f09 
Mean  wind  speed  in  m/s  (solid)  and  spread  in  m/s 
(dashed) . 


Figure  17.  Central  California  Cross  Section: 
Onshore  Flow  fOO. 

Mean  wind  speed  in  m/s  (solid)  and  spread  in  m/s 
(dashed) . 
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Figure  19.  Point  Conception  Cross  Section: 

Along  Coast  Flow  f03. 

Mean  wind  speed  in  m/s  (solid)  and  spread  in  m/s 
(dashed) . 


Figure  22.  Point  Conception  Cross  Section:  12z  f06 
Mean  wind  speed  in  m/s  (solid)  and  spread  in  m/s 
(dashed) . 
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Figure  26.  500mb:  Onshore  Flow  fOO. 

Mean  wind  speed  in  m/s  (solid)  and  spread  in  m/s 
(dashed) . 
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